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Abstract We present results of X-ray spectroscopy for TWA 5, a member of the young TW Hydrae association, observed 
with XMM-Newton. TWA 5 is a multiple system which shows Ha emission, a signature typical of classical T Tauri stars, but 
no infrared excess. From the analysis of the RGS and EPIC spectra, we have derived the emission measure distribution vs. 
temperature of the X-ray emitting plasma, its abundances, and the electron density. The characteristic temperature and density 
of the plasma suggest a corona similar to that of weak-line T Tauri stars and active late-type main sequence stars. TWA 5 also 
shows a low iron abundance (~ 0. 1 times the solar photospheric one) and a pattern of increasing abundances for elements with 
increasing first ionization potential reminiscent of the inverse FIP effect observed in highly active stars. The especially high 
ratio Ne/Fe ~ 10 is similar to that of the classical T Tauri star TW Hya, where the accreting material has been held responsible 
for the X-ray emission. We discuss the possible role of an accretion process in this scenario. Since all T Tauri stars in the TW 
Hydrae association studied so far have very high Ne/Fe ratios, we also propose that environmental conditions may cause this 
effect. 
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1. Introduction 

T Tauri stars are young late-type stars with an age of a few 
Myr, contracting toward the z ero age main sequence phase 
jFeigelson & M ontmerle Hl999t and references therein). They 
are classified in two groups: classical T Tauri stars (CTTSs) and 
weak-line T Tauri stars (WTTSs). This classification is based 
on Ha emission. CTTSs show strong Ha emission (EW > 
10 A). They are still accreting material from their circumstellar 
disk, and broad and asymmetric Ha emission is a direct evi- 
dence of this process. In WTTSs Ha emission is less strong, 
indicating that the accretion process has ended and the star is 
approaching the main-sequence. In most cases CTTSs are also 
characterized by an infrared excess which marks the presence 
of a circumstellar disk. The infrared excess is usually consid- 
ered a prerequisite for accretion, but it does not imply that ac- 
cretion actually takes place; in fact, some WTTSs also show 
an infrared excess although much fainter than in CTTSs. Since 
coeval CTTSs and WTTSs are often observed in the same star 
forming region, the duration of the accretion phase appears to 
be different from star to star. 



One of the signatures of stellar youth is a high X-ray 
emission level. Many star forming regions have been under 
investigation in order to infer the properties of X-ray emis- 
sion from pre-main-sequence (PMS) stars. One of the debated 
questions is whether and how the X-ray emission of accret- 
ing CTTSs and non accreting WTTSs differs. It is conceiv- 
able that the occurrence of the accretion process in CTTSs 
might play a role in determining the different X-ray emission 
characteristics. In fact, the circumstellar disk is thought to af- 
fect the geometry o f the stellar magnetosphere dKoniglll 1 99 lb 
Bouvier et al. 2003). Moreover accreting material may provide 
an alternative heating mechanism for the emitting plasma, al- 
though shock heated plasma cannot attain temperatures higher 
than a few MK. The picture emerging from the analysis of 
low resolution X-ray spectra of PMS stars is that the X-ray 
luminosity of CTTSs is lower than that of WTTSs, and the 
X-ray sp ectra produced by CTT S s appear harder than WTTS 



Tsuiimoto et al. 


2002; 


Flaccomio et al. 


2003; 


IStassun et alJ 


2004; 


Ozawa et 


aj]l2005l). The harder X-ray spectra of CTTSs 
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may be explained with the presence of plasma hotter (T ~ 10- 
100 MK) than that of WTTSs (T ~ 10 MK, Tsuii moto et alJ 
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Table 1. Characteristics of the sample of T Tauri stars, sorted by Ha emission, for which high resolution X-ray spectra have been 
analyzed. In case of resolved multiple systems the stellar mass, spectral type, and bolometric luminosity refer to the component 
responsible for the X-ray emission. 



Name 


Mass 


Spectral 


EW(Hq-) 11 




log(4/L bol ) 


Nl 


References' 1 




(M ) 


Type 


(A) 


(ergs" 1 ) 




(cm- 3 ) 




TW Hya 


~ 0.7 


K7 


-220.0 


1.3 x 10 3U 


-2.7 


~10 13 


1,2, 3,4 


TWA 5 


~ 0.5 


M1.5 


-13.4 


6.7 x 10 29 


-3.1 


< 10" 


2,5,6 


HD 98800 


~ 1.1 


K5 


0.0 


4.1 x 10 29 


-3.8 


< 10" 


2, 7, 8 


PZTel 


~ 1.1 


K0 


0.1 


2.2 x 10 30 


-3.2 


< 10 12 


9, 10, 11 


HD 283572 


~ 1.8 


G5 


1.1 


7.8 x 10 30 


-3.1 




12, 13, 14 



a Negative values of Ho- equivalent width mark an emission line. b X-ray luminosity evaluated in the 6 - 20 A band, using the XMM/MOS 
or Cha ndrafHETGS best fit models pres ente d in the relevant papers. c Densitie s estimated from the O vn and Ne i x triplets. d Data f rom: 
mlBatalha etalJ <2002l) : (2)lReid 120031): ( 3) IStelzer & Schmitil |2004h : (4) lKasmer et"ai]fe002l): (5) this work; (6)ljensen et alj <1998h : (7) 
[ Kastner et alJ fe004l): (8) | Prato et alj feOOll): (9) ICutispoto et all i2002l): (lO)lThatcher & Robinsonl ll993h : (1 D lArgiroffi et alj ^2004^ (12) 
Istrassmeier & Ricd Il998f) : (13) lFemandez & Mirandd ll998h : 7l4) IScelsi et alJ<2005h . " 



120021) . If this is the case, the shock heating mechanism cannot 
be responsible for the X-ray emission in CTTSs. However, it is 
also possible that circumstellar material absorbs the softest part 
of the X-ray ra diation, s imulating therefore a higher tempera- 
ture in CTTSs dStassun et all20 04). 



High resolution X-ray spectra, such as those obtained to- 
day with grating spectrometers on board XMM-Newton and 
Chandra, offer the possibility to reconstruct the emission mea- 
sure distribution (EMD) of the emitting plasma, to measure its 
abundances, and to constrain the electron density N e . These di- 
agnostics help to improve our understanding of the X-ray emis- 
sion from accreting and non accreting young stars. However, 
to achieve a good S IN ratio in these spectra, bright and 
nearby sources ar e needed. The TW Hydrae association (TWA, 
IZuckerman et al-lboOlL and references therein) is one of the 
nearest (~ 55 pc) and youngest (~ lOMyr) star forming re- 
gions and therefore its members are ideal targets for the analy- 
sis of X-ray emission from PMS stars by means of high resolu- 
tion spectroscopy. In the present paper we report on the XMM- 
Newton observation of TWA 5 (CD -33°7795). High resolution 
X-ray spectra of PMS stars have been analyzed in sufficient 
detail so far for only four other stars: TW Hydrae (TW Hya 
or TWA 1), HD 98800 (TWA 4), PZ Tel and HD 283572. 
TW Hydrae and HD 98800, together with TWA 5, are members 
of the TWA; PZ Tel and HD 283572 belong to the /3-Pictoris 
moving group and to the Taurus-Auriga star forming region, 
respectively. 



This paper is organized as follows: in Sect.[2]the principal 
characteristics of TWA 5 and of the other PMS stars used for 
comparison are reported; Sect.[3]presents the main information 
about the XMM-Newton observation of TWA 5 and the methods 
adopted for the data analysis; in Sect. |4] we report the results 
derived, which are discussed and compared with properties of 
other PMS stars in Sect. [5] we draw our conclusions in Sect.|6] 



2. Star Sample 

TWA 5 is a quadruple system located ~ 55 pc from the 
Sun 1 . The primary, TWA 5A, is a triple system: a (X'06 
binary resolved by a daptive optics jMacintosh et alJ l200lt 
iBrandeker et alJ Eo03l) . one of the visual compo nents being 
itself a spectroscopic binary jTorres et alJ 1200 lb . All three 
of them have similar spectral types (Ml. 5). The secondary, 
TW A 5B, is a brown dwarf separated by 2" from the primary 
jLowrance et al]ll999t Iwebb et alJll999h . TWA 5A does not 
show any infrared excess indicat ing no significant amount of 
circumstellar mater ial ( Met chev et al J 120041 [Weinberger et al. 
l2004l lUchida et alJl2004 . On the other hand. iMohantv et all 
( 2003) measured Ha emission typical of accreting PMS stars 
and signatures of outflows, and concluded that at least one of 
the components in the TWA 5 A system is a CTTS. It remains 
currently unclear how the accretion signatures can be recon- 
ciled with the lack of evidence of a disk. Moreover it is un- 
known whether the X-ray emitting component of TWA 5 A co- 
incides with the accreting one. 

Table ^ summarizes the relevant stellar parameters for 
TWA 5 and the other stars, that we will use for comparison 
in our study. 

TW Hya is a single CTT S with enhanced Ha em i ssion 
(equiv alent width ~ 200 A, lAlencar & Batalhal 120021 [Reid 
2003) and strong infrared excess dUchida et alJl2004. Its X- 
ray e mission, observed with Chandra /HETGS ( Kas tner et alJ 
2002) and XMM-Newton dStelzer & Schmittll2004 . shows pe- 
culiar features: the emitting plasma has quite a low tempera- 
ture (log T(K) ~ 6.5), the Ne/Fe abundance ratio is as high as 
a factor 10 in solar photospheric units, and the electron den- 
sity N e , derived from the He-like triplets of Ovu and Neix, is 
~ 10 13 cirT 3 , more than two orders of magnitude above that 
of typical stellar coronae. Based on these peculiarities it was 
suggested that X-ray emission from TW Hya is produced in 
an acc r etion s hock rather than in a corona. On the other hand, 
iDrakd J2005I) has pointed out that the He-like emission line 

1 Only four TWA members have measured Hipparcos distances, 
whose average value, 55 pc, has been assumed as the distance of 
TWA 5. 
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Instrument 


Science 


Filter 


Start 


Exposure 


Count Rate 




Mode 




(UT) 


(ks) 


(cts s' 1 ) 


PN 


Full Frame 


Medium 


2003 Jan 9 03:28:56 


27.9 


1.68 


MOS1 


Full Frame 


Medium 


2003 Jan 9 03:06:55 


29.5 


0.45 


MOS2 


Full Frame 


Medium 


2003 Jan 9 03:06:55 


29.5 


0.46 


RGS1 


Spectroscopy 




2003 Jan 9 03:06:03 


29.7 


0.06 


RGS2 


Spectroscopy 




2003 Jan 9 03:06:03 


29.7 


0.08 



triplets may be affected by photoexcitation due to the UV ra- 
diation field. If this were the case the triplet f/i ratio would 
overestimate the density in the emitting region. 

HD 98800 is a WTTS quadruple system, composed by 
two visual components HD 98800A and HD 98800B (whose 
separation is 0'.'8), each of which is a spec troscopic binary. I t 
has been observed with ChandrafHETGS jKastner et al.l20 04): 
from this observation it emerged that its X-ray emission is due 
mainly to HD 98800A, and it is produced by plasma at tem- 
peratures in the range log 7\K) ~ 6.4 - 7.0, with Ne/Fe ~ 5, 
and low electron density (iV e < 10 I2 crrT 3 ), typical of stellar 
coronae. 

PZ Tel and HD 283572 are two single WTTSs. The X- 
ray spectrum of P Z Tel, gathered with ChandrafHETGS, has 
been analyzed by Argir offi et al.l (l2004). The X-ray spectrum of 
HD 283572, observed with bot h ChandrafHETGS an d XMM- 
Newton, ha s been studied by Aud ard et al.l J2005I) and by 
IScelsi et alJ d2005l) . For both PZ Tel and HD 283572 a typi- 
cal coronal plasma emerged, with temperatures of ~ 10MK, 
and with Ne/Fe ~ 2 - 3 times the solar photospheric ratio. 

3. Observation and Data Analysis 

TWA 5 was observed for ~ 30 ks with XMM-Newton on 2003 
January 9. In Table |2] we report the observation log for all the 
instruments. 

Both EPIC and RGS data have been processed with SAS 
V5.4.1 standard tools. We have extracted EPIC source events 
from a circle with radius of 60", centered on the target posi- 
tion. This extraction circle includes 90% of the source encircled 
energy. Background events have been extracted from an annu- 
lar region around the target with inner and outer radii of 60" 
and 90". The observation is not affected by significant back- 
ground contamination due to solar flares and therefore no time 
screening was required. We have also verified that EPIC spec- 
tra are not affected by significant pile-up. The spectral analy- 
sis has been performed by adopting the Astrophysica l Plasma 
Emission Database (APED VI. 3, ISmith et alJl200lh. which 
assum es ionization equilibrium according to iMazzotta et ail 
Jl99«h . 

3.1. EPIC Data Analysis 

The PN light curve of TWA 5 (Fig. ^ does not show strong 
flare-like events, but an unbinned Kolmogorov-Smirnov test 
applied to the PN photon arrival times yields a probability of 
2x 10~ 7 related to the hypothesis of constant emission. This re- 
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Figure 1. XMM-Newton/PN light curve of TWA 5 with bin size 
of 500 s. 



suit indicates the presence of significant small amplitude vari- 
ability. The observed PN and MOS spectra are shown in Fig. [2 
We have fitted separately the PN and MOS spectra in the en- 
ergy range 0.3 - 7.9 keV, assuming an absorbed optically-thin 
plasma model with three thermal components. We have also 
left as free parameters the abundance of those elements (O, 
Ne, Mg, Si, S, Fe, and Ni) which significantly improved the 
fit, while the abundances of the remaining elements (C, N, Al, 
Ar, and Ca) were tied to the Fe abundance. We have performed 
the fitting by using XSPEC VI 1.3.0. The uncertainty on each 
best-fit parameter, at the 68% confidence level, has been com- 
puted by exploring the x 2 variation while varying simultane- 
ously all the other free parameters. From the PN best-fit model 
we have derived an estimate for the hydrogen column density, 
Nn ~ 3 x 10 20 cm -2 , and the same value has been found as an 
upper limit from the analysis of the MOS spectra. This result 
indicates that the spectra of TWA 5 do not suffer strong absorp- 
tion. The der i ved A^ h value is compatible with that assumed by 
jjensen et alJ dl998t> . which agrees with the negligible extinc- 
tion toward the TWA region. The results obtained from the PN 
and MOS spectral fitting are reported in Table [5] 

3.2. RGS Data Analysis 

The RGS1 and RGS2 spectra of TWA 5 are shown 
in F ig. 121 The analysis has been perf ormed using ISIS 
faouck & Denicolal2000h and PINTof ALE JKashvap & Drake! 
2000). Our approach is to derive EMD and abundances start- 
ing from the line flux measurements. It is known that the line 
spread function of the RGS spectra is characterized by large 
wings, making it difficult to identify correctly the continuum 
level and therefore to measure line fluxes. In order to obtain 
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energy (keV) energy (keV) 



Figure 2. EPIC spectra (PN in upper panel, MOS 1 and MOS2 in lower panel) of TWA 5 with best-fit 3-T models superimposed. 
The lower section of each panel contains residuals. Best-fit parameters are reported in Table[3] 




5 10 15 20 25 30 35 

Wavelength (A) 



Figure 3. Observed and predicted RGS spectra (rebinned with a wavelength bin of 0.03 A) of TWA 5. 



accurate line flux measurements we have evaluated the contin- 
uum level by performing a global fit of the RGS 1 and RGS2 
spectra. We have adopted a model composed of three isother- 
mal components with variable abundances of C, N, O, Ne, Mg 
and Fe. The continuum predicted on the basis of this best-fit 
model has been used to measure the fluxes of the strongest 
RGS emission lines. To improve the spectral S [N ratio we have 
measured the line fluxes by fitting simultaneously RGS1 and 



RGS2 spectra rebinned with a 0.03 A wavelength bin. These 
line fluxes are reported in Table|5f ■ 

We have reconstructed the EMD and element abundances 
with the Markov-Chain Monte Carlo (MCMC) method of 

2 Table [5] is available at the CDS and it contains the following in- 
formation for each observed line: observed and predicted line wave- 
length (Cols. 2 and 3), element and ionization state (Col. 4), electronic 
configurations of the atomic levels (Col. 5), temperature of maximum 
emissivity (Col. 6), observed line flux (Col. 7). 
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iKashvap & Drakel dl998h applied to the measured line fluxes. 
This method performs a search in the EMD and abundances 
parameter space with the aim of maximizing the probability of 
obtaining the best match between observed and predicted line 
fluxes. Some of the main advantages of this method are that 
it does not need to assume a particular analytical function for 
the EMD, and that it allows to estimate uncertainties on each 
EMD and abundance value. On the basis of the formation tem- 
perature of the selected set of lines, we have adopted a temper- 
ature grid ranging from log T(K) = 6.0 to log T(K) = 7.3, with 
resolution A log T(K) = 0.1, over which to perform the EMD 
reconstruction. We have assumed a hydrogen column density 
A^h = 10 20 crrT 2 , compatible with the values derived from the 
analysis of the EPIC spectra (see Sect. 13. 1> . In Fig. @] we show 
the comparison between the observed line fluxes and those pre- 
dicted on the basis of the EMD and abundances derived with 
the MCMC method. 

Since line fluxes depend on the product of the EMD with 
the element abundances, the adopted method provides the 
EMD scaled by the Fe abundance, and the abundance ratio of 
each element with respect to Fe. However, it is worth noting 
that the continuum emission, depends strongly on the amount 
of emission measure, and weakly on the absolute abundances 
of elements heavier than He. In fact the continuum emission is 
due to three processes: bremsstrahlung radiation, radiative re- 
combination and two-photons emission. For the temperatures 



Table 3. TWA 5 best fit paramenters. 



PN 



MOS 
EE" 



RGS 



Abundances 1 (A x /A x q) 



c 


= Fe 


= Fe 


26+ 030 

"•^"-0.08 


N 


= Fe 


= Fe 


55 +0 - 61 

U - J -0.18 


O 


15 +on 

"■ iJ -0.07 




31 +0.29 


Nc 


34 +0 - 27 

"•-"-0.16 


0.93*$ 


92 +0 - 74 


Mg 


07 +015 

"•"'-0.06 


14 +020 

ul -0.13 


51+° 16 

-0.38 


Si 


021+0.20 


19 +022 


= Fe 


S 


32 +036 


36 +0 38 

-0.32 


= Fe 


Fc 


05 +0 04 

"• UJ -0.02 


09 +0 06 

-0.05 


0.1 


Ni 


1.08+J;' 7 , 


06 +l 31 


= Fe 



Temperature b (K) 


log T, 
logT 2 
log T 3 


6 17 +0 - 56 

'-0.20 

6 70 +017 

7 97+O.09 
'•^'-0.06 


6 56 +0 14 

u ' 7 "-025 

7 97+0.33 
' 1 -0.09 


Emission Measure b N e NyiV (cm 3 ) 


log EM i 
log EM 2 
log EM 3 


52 71+0-™ 

' x -0.63 

53 21 +0 - 26 

JJ -0.31 
S9 79+0.10 


S? 4.S+0- 60 

J J -0.37 

52 76 +0 34 

JZ, -' U -0.70 

52 64 +015 

JZ "" -0.30 


Column Density b (10 20 cirT 2 ) 


N H 


2.8! 3 ;^ 


< 3.3 =1 


Best Fit Statistics 


^red 

d.o.f. 


0.89 
404 


0.99 
285 


p (x 2 >xtJ 


94% 


54% 



a Solar photospheric abundances are from lAnders & Grevessel 
<1989h . 

b All the uncertainties correspond to the 68% confidence level. 



4 - a c • n o o 

2 7 

o r « «■ 

-2 r 

-4 1 



* Ne + Mg o Fe 



6.0 



6.2 



6.4 6.6 6.8 
log T (K) 



7.0 



7.2 



Figure 4. Comparison between observed and predicted line 
fluxes in the RGS spectra. The predicted values are obtained us- 
ing the models derived with the MCMC method (see Sect. 13. 2\ . 

involved in the plasma of TWA 5 the main contribution to the 
continuum is due to bremsstrahlung radiation which depends 
very weakly on the heavy element abundances. Therefore, af- 
ter performing the MCMC reconstruction, we have considered 
several models assuming different absolute Fe abundances and 
therefore different global scaling factors for the EMD distri- 
bution. For each of these models we have compared the pre- 
dicted and observed continuum levels. Since it is hard to iden- 
tify correctly the continuum level in RGS spectra, as already 
mentioned above, we have also compared the observed and 
predicted total emission (spectral lines + continuum) as a fur- 
ther check. In Table 0] we report the explored Fe abundances, 
and the corresponding total counts for the simulated spectra 
A'pred, to be compared with the observed total number of counts, 
A^obs- With this procedure we have determined the absolute Fe 
abundance, and therefore the absolute position of the EMD and 
the absolute abundances of all the other elements. The result- 
ing Fe abunda nce is 0.1 tim es the solar photospheric value of 
lAndersfe Grevessel lll989h . with an uncertainty smaller than 
a factor 2. As a final cross-check, we have verified that the 
predicted continuum level agrees with the continuum used for 
the line flux measurements. The abundances resulting from the 
RGS analysis are reported in Table|5] 

4. Results 

The X-ray luminosity of TWA 5, computed in the interval 
6 - 20 A from the best-fit models of PN, MOS and RGS spec- 
tra, is 8.3, 6.7 and 6.7 x 10 29 erg s _1 , respectively. The derived 
luminosities are compatible wi thin th e best-fit parameter er- 
rors. As shown bv lTsuboi etaTl lExB). from a Chandra/ ACIS- 



Table 4. Total counts (lines + continuum) in the interval 10- 
30 A band vs. Fe abundance. 















/Vobs 


Fe/Fe a 


0.05 


0.07 


0.10 


0.15 


0.20 




RGS1 


2298 


1841 


1498 


1232 


1099 


1377 


RGS 2 


2886 


2379 


1998 


1702 


1554 


1874 



a Abundance referred to the solar photospheric value from 
lAnders & Grevessd JT989I) . 
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Ni Mg Fe Si S C N 
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Figure 5. Emission measure distribution of TWA 5 obtained 
from the RGS line fluxes with the MCMC method, along with 
the values obtained with the 3-T best-fit model of PN and MOS 
spectra. The uncertainties correspond to 68% statistical confi- 
dence level. 

S observation, the X-ray emission is essentially due to the pri- 
mary TWA 5A. In fact the Chandra observation was able to 
resol ve the brow n dwarf TWA 5B from the primary TWA 5A, 
and lTsuboi et al J measured for TWA 5B an X-ray luminosity of 
4x lO^ergs" 1 . 

4.1. Emission Measure Distribution 

In Fig.|5]we report the EMD vs. temperature derived from the 
EPIC and RGS spectral analysis of TWA 5 3 . All the instru- 
ments detect the strongest thermal component at log T(K) ~ 
6.7 - 7.0, but the RGS spectra are not able to probe the hottest 
plasma component at log7\K) ~ 7.3, detected by EPIC. The 
reason of this result is the different effective area of EPIC and 
RGS in the hardest part of the X-ray spectra (E S 2keV). In 
principle, the high-temperature tail could be probed by exploit- 
ing a number of Fe xxii-xxiv lines, which fall in the wavelength 
region 10.7 - 11.8 A covered by RGS, but the emissivity of 
these lines is relatively low and the RGS resolution too poor 
for this purpose. 

4.2. Abundances 

In Fig. |6] (and in Table we show the element abundances, 
in solar photospheric units dAnders & Grevessell989l) . derived 
from the spectra obtained with each instrument. The elements 
are sorted along the abscissa by increasing values of first ion- 
ization potential (FIP). The abundances of C and N, which have 
their strongest emission lines in the low-energy part of the ob- 
served spectral range (A ~ 25-35 A, or E ~ 0.35 -0.5 keV), are 

3 We note that the MCMC method explores preferentially EMD 
bins which are best constrained by the selected emission lines. Since 
error estimation depends on the quality of the sampling, statistical 
uncert ainties are estim ated only for those EMD bins explored many 
times l Kashvap & Dra kell998h . 



Figure 6. Abun dances of TWA 5, with re spect to solar photo- 
spheric values dAnders & Grevessel[l9 89). obtained using the 
RGS line fluxes with the MCMC method, along with the values 
obtained with the 3-T best-fit model of PN and MOS spectra. 
The uncertainties correspond to the 68% confidence level, ex- 
cept the error bar of Fe derived from RGS data, which has been 
obtained with a different procedure, as discussed in Sects. l3~2l 
andl4~2l 

derived only from the RGS. On the other hand, the abundances 
of Si and S are estimated only from the EPIC spectra since their 
H-like and He-like lines fall at high energies, and they cannot 
be constrained by the RGS. Note that the Fe abundance derived 
from RGS data has been estimated with a procedure (Sect. 13. 2l 
which does not allow to obtain a formal statistical uncertainty, 
however we are confident that it cannot be off by more than a 
factor 2, as explained in Sect. 13.21 We note that the abundances 
derived from different instruments are compatible within statis- 
tical uncertainties. However some systematic differences come 
out, as briefly discussed in Sect. 14.41 From the derived abun- 
dances it emerges that the X-ray emitting plasma of TWA 5 is 
metal depleted. 

4.3. Electron Density 

We have evaluated the plasma electron density from the analy- 
sis of the O vii He-like triplet. The other He-like triplets which 
fall in the RGS spectral range were either too weak (Nvi, 
Mg xi, Si xm), or heavily blended with other strong lines (Ne ix) 
to be analyzed. In Fig.[7]we show the RGS 1 spectrum (rebinned 
with a 0.03 A bin size) in the O vn triplet region with superim- 
posed the best-fit curves. The measured ratio of forbidden / 
and intercombination i line fluxes is 3.8 +2.5 (see Table [5]for 
line fluxes), which yields an upper limit of ~ 10 11 cm' 3 for N e , 
adopting the predicted // i ratios of Smi th et alJ (|2001 ). 

4.4. Comparison between Different Models 

We have performed separate analyses of the PN, MOS and RGS 
spectra for several reasons. The main reason is that an approach 
based on fluxes of selected individual lines, measurable only in 
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21.4 21.6 21.8 22.0 22.2 
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Figure 7. Ovn triplet in the RGS1 spectrum of TWA 5, re- 
binned with a wavelength bin of 0.03 A, with best-fit line pro- 
files: individual lines (dotted line) and their sum (solid line). 



the RGS data, provide us with the most reliable results for the 
element abundances and for the plasma EMD. Moreover, inde- 
pendent analyses of EPIC data offer the opportunity to compare 
the results of the different XMM-Newton instruments. These 
comparisons allow us to investigate the robustness of each mea- 
surement, and therefore they are useful to test the reliability 
of results based on EPIC data only, in the broader contest of 
observations of X-ray coronal sources with no high resolution 
spectrum available. 

As already mentioned in Sect. 14.21 abundance estimation 
from different XMM-Newton instruments turns out to be quite 
robust, at least within the statistical uncertainties of a typi- 
cal XMM-Newton observation (~ 30 ks exposure time, in the 
present case). However, the abundances of Mg, Fe, O, and Ne 
obtained from the PN fitting are systematically lower than the 
corresponding values based on the analysis of the MOS and 
RGS spectra, which agree among themselves. On the other 
hand, the emission measure values derived from the PN anal- 
ysis are larger than those obtained from the MOS and RGS 
spectra, and the X-ray luminosity predicted by the PN model is 
~ 20% higher than in the other two cases. Although, all the dif- 
ferences are within the statistical uncertainties, it is conceivable 
that the higher spectral resolution of the MOS detector, with re- 
spect to the EPIC/PN, allows to disentangle better the contribu- 
tions of lines and continuum, and therefore to constrain the ab- 
solute values of EMD and abundances. Moreover, EPIC/MOS 
and RGS share the same X-ray telescopes and hence their 
cross-calibration is better determined, while residual calibra- 
tion problems of the PN instrument may cause the differences 
in the fitting results described above. 

EPIC models are able to provide a good global description 
of the source plasma, but limited to models with few free pa- 
rameters, while RGS spectra allow to derive a more detailed 
model but - with the available signal to noise ratio - they fail 
to detect plasma components with temperatures higher than 
T > 10MK, due to the smaller energy range covered by RGS 
with respect to EPIC. 



In order to cross-check the three models we have compared 
each of them with the spectra of different XMM-Newton in- 
struments and we have computed reduced^ 2 values. The best- 
fit 3-T models of PN and MOS globally describe RGS spec- 
tra reasonably well (y 2 ed = 1.8 and ^ 2 ed = 1.3 respectively, 
with 817 d.o.f.), but not as well as the line-based EMD model 4 
(*red = 10 with 817 d.o.f.). As already noticed the RGS model 
misses the higher temperature components and hence it under- 
estimates the PN and MOS spectra at high energy. Finally, both 
RGS and MOS models show some disagreement with the PN 
spectrum in the low energy range. 



5. Discussion 

In this section we discuss the results obtained for TWA 5, in 
terms of EMD, abundances and density, and compare them 
with those for the CTTS TW Hya, and the other WTTSs in 
our sample (Sect. |2}. We stress that all results are based on 
high resolution X-ray spectra. It must be recalled that TWA 5 A 
is a triple system, and so far we have not been able to deter- 
mine whether the X-ray emission and the accretion signatures 
emerge from the same star. 

The analysis of both the EPIC and RGS data has shown that 
the X-ray emission of TWA 5 is mainly produced by hot plasma 
(T ~ 10 MK). The analysis of the Ovn triplet has indicated a 
typical coronal electron density (N e < 10 11 cm" 3 ). These char- 
acteristics are similar to those found in WTT Ss I Kastner et alJ 
l2004UArgiroffi et alJ20o4lScelsi et al.l2005l) and magnetically 
active late-type main sequence sta rs when high reso lution X- 
ray spectroscopy is used (see e.g. iNess et al.ll2004l) . Peculiar 
features of TWA 5 are its very low metallicity (Fe/FeQ ~ 0.1) 
and its extremely high abundance ratio Ne/Fe ~ 10. Such high 
values for the Ne/Fe have been observed only in a few very 
active stars (HR 1099, UX Ari, II Peg) having average coronal 
temperatures larger than those of TWA 5 (see below), and it has 
been ascribed to the so-called inverse FIP effect (see discussion 
below). On the other hand, a low metallicity (Fe/FeQ ~ 0.2) 
and the same Ne/Fe ratio have been measured for TW Hya, the 
only unambiguous CTTS studied so far at high spectral reso- 
lution in X-rays. Hence, it is an interesting issue why TWA 5 
shares the same chemical peculiarities with TW Hya, in spite 
of having other thermal characteristics. 

We recall that TW Hya presents spectral characteristics 
compatible with a model of X-ray emission driven, or at 
least affected, by the infalling accretion stream JKastner et alJ 
2002| IStelzer & Schmittll2004l) . In fact, all the X-ray proper- 
ties of TW Hya (e.g. its low plasma temperature, high density, 
and metal depletion) suggest that the emitting plasma forms 
in the shock region where the infall streams reach the stel- 
lar surface. Among the peculiar characteristics of TW Hya 
the very low abundances of all the metals in the emitting 
plasma appea r to be co mpatible with the accretion scenario. 
In fact, IStelzer &"S chmitt proposed that Fe and other heavy 
elements in the accretion disk condense into dust grains (see 



4 This is not obvious since our RGS EMD model has been obtained 
from the analysis of selected line fluxes and not from a global spectral 
fitting. 
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Figure 8. Ne/Fe abundance ratio in solar photospheric units 
jAnders & Grevessdll989l) for the PMS stars listed in Table [JJ 
and sorted by Ha equivalent width. We have assumed that the 
abundance uncertainties of HD 98800 are equal to those of 
TW Hya, because the observed Chandra spectra of these two 
stars have similar S /N. 



Figure 9. Ne/Fe abundance ratio in solar photospheric units 
dAnders & Grevessel 19891 for the PMS stars listed in Table[TJ 
vs. the X-ray to bolometric luminosity ratio. We have assumed 
that the abundance uncertainties of HD 98800 are equal to 
those of TW Hya, because the observed Chandra spectra of 
these two stars have similar S IN. 



e.g. Savage & Sembachl l 19961) which possibly settle into the 
disk midplane, while other elements like N, remain in the gas 
phase. Neon and other noble elements, should also not remain 
locked onto dust grain s, but rather be part of the gas phase 
(Frisch&Slavi nll2003h . 



Since the accreting material is largely 
composed by gas rather than dust (Takeuc hi & Linl l2002). the 
accreting stream is expected to display a high Ne/Fe abundance 
ratio. This material falls onto the stellar surface and there, 
heated to temperatures of few MK by the ensuing shock, pro- 
duces X-ray radiation revealing its anomalous chemical com- 
position. The intriguing point is that TWA 5 presents exactly 
the same abundance ratios of TW Hya, and in particular a 
Ne/Fe ~ 10, but lacks all the other indications for accretion- 
related X-ray emission. In Fig. [8]we show the Ne/Fe ratio for 
the PMS stars in our sample. The stars are sorted along the ab- 
scissa by increasing value of Ha equivalent width, in order to 
separate the accreting CTTSs, on the left part of the diagram, 
from the non accreting WTTSs. This plot suggests that CTTSs 
tend to have Ne/Fe higher than WTTSs. 

As already hinted above the Ne/Fe ratio could be influenced 
also by FlP-related effects: in the solar corona, and in partic- 
ular in long-lived active regions, and in late type stars with 
low activity levels, abundances of elements with low FIP ap- 
pear to be enhanced with respect to the high FIP elements (see 
iFeldman & Widingll2003l and references therein), using pho- 
tospheric abundances as a reference. On the other hand, more 
active stars present an overabundance of high FIP elements 
with r espect to low FIP elements, the so called inverse FIP 
effect JBrinkman et aljEool iDrake et alJl200lt lAudard et alJ 
2003). Early models to explain the FIP effect involve ion- 
neutral fractionation in the chromosphere (Geiss 1982; Mever 
Il996l) . but they do not provide a satisfactory explanation for 
the selecti ve enh ancement of some elements in th e corona 
(see)Gudel 20041 for a recent review). Most recently, Laming 



(2004) has proposed a new model which tries to explain both 
a FIP and an inverse FIP effect as a result of the ponder- 
motive forces related to chromospheric Alfven waves acting 
on ions of different species. The Ne/Fe ratio is a good indi- 
cator of the coronal abundance pattern since Ne has a high 
FIP value (21.6 eV), while Fe is a low FIP element (7.9 eV), 
and strong lines from both elements have close wavelengths 
at similar coronal temperatures. Stars with high activity level 
usually show Ne/Fe ~ 1-5, and only few active binaries 
present Ne/Fe ~ 10 JBrintoanetalJl200lUDrakeetalJ|200ll: 
iHuenemoerder et alJl200lt lAudard et alJl2003h . lGudeiT (l2004) 
shows that the Ne/Fe ratio tends to increase for increasing 
average coronal temperature, with the above extreme value 
reached by stars with T c ~ 15 MK. For comparison, TWA 5 
has T c « 9MK and sta rs of comparable temperature in the 
sample studied by Gudel show Ne/Fe in the range 1-5. 

To explore further this inverse FIP effect scenario we have 
plotted the Ne/Fe ratio for the PMS star sample vs. Lx/^boi 
in Fig. [5] In fact active stars do show a c orrelation between 
coronal a bundances and th e activity l evel dSingh et al]li993 
iGiidel et alJl2002t lAudard et al]l2003l) . If the differences in 
Ne/Fe ratio among the stars in our sample were caused by a 
similar FlP-related effect we would expect to see a correlation 
between Ne/Fe and Lx /Lboi- For comparison purposes, we have 
included in the plot also TW Hya, even if its X-ray emission is 
likely not due to coronal activity. This plot does not show any 
clear trend, even if we do not consider TW Hya. This result 
might be due to the small number of PMS stars studied so far 
with high resolution X-ray spectroscopy, and to the fact that 
most of these stars are in the saturated emission regime where 
Lx/^boi ~ 10~ 3 . If we insist that an inverse FIP effect is respon- 
sible for the observed Ne/Fe ratio of TWA 5, it still remains 
unclear why stars with similar characteristics (age, plasma tem- 
perature, Lx/Lbol) do show Ne/Fe values which differ by about 
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a factor 10, as in the case of TWA 5, PZ Tel and HD 283572. 
Hence we argue that TWA 5, and even more clearly TW Hya, 
appear to be outliers with respect to other active stars. 

In conclusion we can tentatively depict three different sce- 
narios in order to interpret the characteristics of the X-ray emit- 
ting plasma in TWA 5. 

Since TWA 5 appears to contain a CTTS ( M ohantv et alJ 

120031) . the high Ne/Fe might be due to an accretion process, 
as already suggested in the case of TW Hya. In this sce- 
nario, the X-ray emission from TWA 5 should be produced 
by shock heated plasma at the base of the accretion column. 
However, shock temperatures are expected to be lower than 
the values derived from the X-ray spectrum of TWA 5, and 
this occurrence is not in favor of accretion-related X-ray emis- 
sion. This scenario is also questioned by the analysis of X- 
ray emission from CTTS s and W TTSs in the L1551 region, 
discussed by Favat a et al 1 il2003h : they derive Ne/Fe ~ 4 for 
the three WTTSs and no indication of high Ne/Fe for the two 
CTTSs in their stellar sample 5 . Most recently the analysis of 
the XMM-Newton/PN spectrum of the CTTS BP Tau revealed 
a hot plasma, while the O vn lines suggested a high electron 
density ( Sch mitt et all200 5). These results on BP Tau indicate 
that shock heated and coronal plasma may be both present in 
CTTSs. 

The second scenario is based on the consideration that 
TWA 5 has \og{Lx I L\, \) ~ -3, at the saturation level for ac- 
tive stars. Therefore, the high Ne/Fe ratio may be related to the 
same mechanism which produces the inverse FIP effect in the 
coronae of other active stars. Under this hypothesis the accre- 
tion process does not play a major role in the X-ray emission 
of TWA 5, which is instead produced by magnetically confined 
hot plasma. However the Ne/Fe ratio of TWA 5 appears to be to 
high by a factor 2-5 with respect to stars with similar average 
coronal temperature. 

Finally, we note that both TW Hya and TWA 5 belong to 
the same young association, and share the same value of Ne/Fe. 
Therefore, the third hypothesis is that their anomalous abun- 
dances originate from the molecular cloud from which the two 
stars formed. Such a scenario, in which the measured abun- 
dances are related to those of the primordial material implies 
that the molecular cloud was Fe depleted. In order to con- 
firm or reject this hypothesis the abundances of other mem- 
bers of the TWA need to be determined. Note that in the case 
of HD 988 00, a member of TW A, the ratio Ne/Fe ~ 5 was 
derived by iKastner et alJ J2004I) from a spectrum affected by 
low S/N ratio which did not allow these authors to perform a 
detailed EMD analysis. As a consequence the derived Ne/Fe 
ratio is uncertain since it depends strongly on the EMD shape. 

6. Conclusions 

We have analyzed the EPIC and RGS data of the CTTS TWA 5 
inferring the emitting plasma characteristics: the X-ray emis- 
sion reveals a hot plasma (T ~ 10 7 K) with low electron density 
(N e < 10" cirT 3 ) and low metallicity (Fe/Fe© ~ 0.1). These 

3 However, the results on L155 1 region are based on low resolution 
EPIC spectra, and therefore they may not be directly compared to ours. 



findings suggest that X-rays may be generated by magnetically- 
confined coronal plasma strongly influenced by an inverse FIP 
effect. However stars with coronal temperatures comparable 
with that of TWA 5 show lower Ne/Fe ratios (Gudel 2004). 
The Ne/Fe ~ 10 abundance ratio measured for TWA 5 leaves 
open the issue of the X-ray production mechanism, since the 
same Ne/Fe has been measured for the CTTS TW Hya, where 
this result has been interpreted as evidence for the shock heated 
accreting material as responsible for the X-ray emission. An al- 
ternative explanation we propose is that the peculiar abundance 
ratio could be a characteristics of the primeval gas from which 
all members of the TWA formed. 
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8r 
oc 




1 ^ 74.^8 


Fe 


XIX 


9 „29 3/2 ry,a j 1 p _. tJtJ In 
zs Lp \ u)ju r 3 * zs z^y i>/2 










Qa 




14. 9080 


Fe 


XVIII 


is zs z/^i/2Z/^3^2 J "5/2 ^ zs z/^ -T3/2 


6 QO 


1Z..Z. 




jL.XJ 






1 /I Of\Qf\ 
14.ZUSU 


Fe 


XVIII 


O n2o ^4/ 1 r»\'2^7 2 ft . o n2o ^.5 2n 

Zs Z/? ( L» 5 /2 — > ZS Zp ^3/2 


a Qr\ 








yc 




14.ZJ0U 


Fe 


XVIII 


1 P 2o r 2o „ o ^3 o^j . /-> „2o ^,5 2p 

IS Zs Lp\i2Lp il2 idsi2 -* LS Lp Pi/2 


a Qr\ 








10a 


1/1 CO 

14.58 


14.4856 


Fe 


XVIII 


o 2o 4/1 r\\ o j 2c . o 2o 5 2n 

2s 2p ( D)3a S 1/2 — > 2s 2/? P1/2 


/■ r\r\ 

6.90 


1 A C 

10.5 


± 


2.5 


1UD 




i /i cnc/r 
14.JUJ0 


Fe 


XVIII 


1 ^20^20^.2 «2 j . 0^20^.5 2d 
IS ZS Lp l/2 Lp 3/2 id3/2 -* ZS Z/? ^3/2 


O.oU 








10c 




14.5340 


Fe 


XVIII 


2o 4/3 TVii j 2r 9 5 2ti 

2s / -2p*{ 3 P)3d F 5 /2 — > 2s' L 2p J P3/2 


6.90 








1 a j 

10d 




14.5710 


Fe 


XVIII 


2s / -2p*( 3 P)3d P3/2 — > 2s^2p J P3/2 


6.90 








ll 


15.01 


15.0140 


Fe 


XVII 


--\ 'J r\ S / 9 / 1\ t 1 n /-\ 9 /^ A I n 

2s z 2p D ( £ P)3d l P\ —> 2s z 2p OL So 


6.70 


19.5 


± 


3.3 


12a 


15.17 


15.1760 


O 


VIII 


A 1 rt 1 2 c 

4p Z /P 3 /2 -» 1« S 1/2 


6.50 


7.0 


+ 


2.5 


12b 




15.1765 


O 


VIII 


4/? 2 Pi/ 2 -» ls 2 Si /2 


6.50 








1ZC 




i c i qqa 


Fe 


XIX 


1 „20 nO r.2 O r.2 1 r, v O P 9 ,,5 3 D 

IS LSLp^ijLp^.jpS — > LSLp r2 


o.yu 








13 


15.28 


15.2610 


Fe 


XVII 


2s / -2p 2 { A P)3d D\ -> 2s z 2p° 'So 


6.70 


3.7 


+ 


2.3 


14a 


15.96 


1 f f\f\ A f\ 

16.0040 


Fe 


XVIII 


2s z 2/? 4 ( J P)3s 2 P 3/2 -> 2s 2 2/? 3 z P 3/2 


6.80 


12.0 


+ 


3.9 


14b 




16.0055 


O 


VIII 


3/? 2 P 3 /2 -» F? 2 Si /2 


6.50 








14c 




16.0067 


O 


VIII 


O 2 75 1 2 C 

3? Pl/2 -» IS Sl/2 


6.50 








15a 


1 a ao 

16.02 


16.0710 


Fe 


XVIII 


/-> 2 ^ 4/3 7"*\ o 4 n a ?a S 2n 

2s z 2/? 4 ( J P)3s 4 P 5/2 -> 2s z 2/5 3 z P 3/2 


6.80 


19.8 


+ 


4.3 


15b 




16.1100 


Fe 


XIX 


1 9^ 9^ 9 o ^> S ^ / 1 

ls z 2s z 2^i /2 2^^ /2 3/?i/2 2s2/? 3J P 2 


6.90 








16 


16.74 


16.7800 


Fe 


XVII 


2s 2 2p 5 ( 2 P)3s 'Pi -> 2s 2 2/? 6 'So 


6.70 


13.3 


+ 


2.8 


17a 


17.06 


17.0510 


Fe 


XVII 


2s 2 2p 5 ( 2 P)3s 3 Pi -> 2s 2 2p 61 S 


6.70 


22.5 


+ 


3.6 


17b 




17.0960 


Fe 


XVII 


2s 2 2p 5 ( 2 P)3s 3 P 2 -* 2s 2 2p 61 S 


6.70 








18 


17.58 


17.6230 


Fe 


XVIII 


2s 2 2p 4 3p 2 P 3 /2 -> 2s2/? 62 Si /2 


6.80 


6.7 


+ 


2.4 


19 


18.66 


18.6270 


O 


VII 


ls3p 'Pi Is 2 'S 


6.30 


5.7 


± 


2.8 


20a 


1 o a/c 

18.96 


18.9671 


O 


VIII 


o 2 n 1 2 c 

2/? z P 3/2 -» ls z Si /2 


6.50 


1 19.7 


± 


6.8 


20b 




18.9725 


O 


VIII 


2/? 2 P 1/2 -» ls 2 Si /2 


6.50 








21 


21.60 


21.6015 


O 


VII 


ls2p 'Pi -> Is 2 'So 


6.30 


39.8 


± 


6.9 


22 


21.80 


21.8036 


O 


VII 


ls2/? 3 Pi -> Is 2 'So 


6.30 


6.4 


± 


3.9 


23 


22.08 


22.0977 


O 


VII 


ls2s 3 S! ^ Is 2 'So 


6.30 


24.5 


+ 


5.5 


24a 


24.79 


24.7792 


N 


VII 


2p 2 P V2 -» ls 2 Si /2 


6.30 


20.8 


+ 


4.3 


24b 




24.7846 


N 


VII 


2/? 2 Pi/ 2 -» ls 2 Si /2 


6.30 









12 C. Argiroffi et al.: XMM-Newton spectroscopy of TWA 5 



Table 5. continued. 



Label 


^obs 


^pred 


Ion 




Transition 




(F 


± 


o-F) c 




(A) 


(A) 






(upper — > lower) 


(K) 




(io- 


6 ph s _1 cm -2 ) 


25a 


33.73 


33.7342 


C vi 


2p *P 3 /2 


-» 1*^1/2 


6.10 


22.5 


± 


4.8 


25b 




33.7396 


C vi 


2 P 2 Pi/ 2 




6.10 









a Observed and predicted (APED database) wavelengths. In the cases of unresolved blends, identified by the same label 

number, we list the main components in order of increasing predicted wavelength. 
b Temperature of maximum emissivity. 

c Line fluxes with uncertainties at the 68% confidence level obtained by fitting simultaneously RGS1 and RGS2 spectra. In 
the cases of unresolved blends, identified by the same label number, we report only the total flux of the blended lines. 



